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We present a simple exact formula for the three point velocity correlation functions in two dimensional 
turbulence which is valid on all scales and which interpolates between the direct and inverse cascade regimes. 

f^ . As expected, these correlation functions are universal in these extreme regimes. We also discuss the effect 

Cn ' of anisotropy and friction. 
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The aim of this short note is to write down an exphcit formula for the three point velocity correlation 
fmictions in 2d turbulence. See eqs.(p|) below and its consequences. This formula differs from the usual 
Kolmogorov formula, cf eg ref. |l| , by the fact that it incorporates the existence of two inertial ranges which 
respectively correspond to the inverse energy cascade and the direct enstrophy cascade. Although expected 
on scaling ground, this formula, and its simple proof, was surprisingly never spelled out in the turbulent 
literature. It is nevertheless one of the rare exact results on 2d turbulence. We thus feel that it was worth 
making it more public. 

Kraichnan's scaling theory. A special feature which distinguishes two dimensional from three dimen- 
sional fluid mechanics is the conservation of vorticity moments in the inviscid limit. As first pointed out by 
Kraichnan in a remarkable paper [gj , this opens the possibility for quite different scenario for the behavior of 
turbulent flows in two and three dimensions. In two dimensions the inviscid Navier-Stokes equation admits 
two quadratic conserved quantities: the energy J ^ and the enstrophy J ^ with u the velocity and w the 
vorticity. As argued by Kraichnan, if energy and enstrophy density are injected at a scale Li, with respective 
rate e and e^ ~ ^L~^, the turbulent system should react such that the energy flows toward the large scales 
and the enstrophy towards the small scales. As this energy flow is quite the opposite to the one involved 
in Kolmogorov's picture for 3d turbulence, one usual refers to the infrared energy flow as the inverse cas- 
cade and to the ultraviolet enstrophy flow as the direct cascade. The fact that energy has to escape to the 
large scales may be understood from the fact that in absence of forcing the time variation of the energy is 
dt J \ = —1^ J ^ with i^ the viscocity. It thus vanishes in the inviscid limit i^ ^ if the enstrophy remains 
finite and the energy cannot be dissipated at small scales. 

In the (IR) inverse cascade, scaling arguments lead to Kolmogorov's spectrum, with E{k) '^ e^/^^ k'^/^ 
for the energy and (AM)(r) '^ {erY'^ for the variation of the velocity on scale r. 

In the (UV) direct cascade, scaling arguments give Kraichnan's spectrum with E{k) ^ ej fc~^ for the 
energy and {/\u){r) '-^ (e^r'^)^/'^ for the velocity variation. 

Of course the direct and inverse cascade have been extensively analysed, both numerically, see eg. [0 
and refs. therein for an (incomplete) sample of references, and theoretically, see eg. 0, p[ and refs therein 
for a few relevent references, some of which discussing logarithmic corrections to Kraichnan's scaling. More 
recently, the inverse cascade has been observed experimentally, as described in ref. H . Within experimental 
precision it shows no deviation from Kolmogorov's scaling. 

A model and its hypothesis. As usual, to statistically model turbulent flows we consider the Navier- 
Stokes equation with an extra forcing term. Let u^{x,t) be the velocity field for an incompressible fluid, 
V • w = 0. In two dimensions the incompressibility implies that u{x, t) derives from a stream function <I> such 
that u'' — ekjdj^ with ekj the antisymmetric tensor. The Navier-Stockes equation reads: 

dtu^ + {u ■ V)u^ - lyV^u^ = -V^p + P (1) 

with p the pressure and /(x, t) the external force such that V • / = 0. We choose the force to be Gaussian, 
white-noice in time, with zero mean and two point function: 

{f[x,t)f\y,s)) = a\x~y)5[t~s) (2) 

where C^^{x), with V-' C^*'{x) = 0, is a smooth function varying on a scale Li, fastly decreasing at infinity 
and regular at the origin. The scale Li represents the injection length. We shall assume translation, rotation 
and parity invariance, unless otherwise specified. Let C{x) = trC{x). Its Taylor expansion at the origin 
will be denoted as: C{x) = 2e — e^, r^/2 + • • ■ with r^ — x^xt- The transversality condition V^ C^'^{x) = 
ensures that C = V^Q'^ with 0*'(a;) — Jx^ — e^ x'^r'^ /8 -I- • • • at short distances. Physical interpretation of e 
and Cu, will be given later. 

The vorticity w, with oj — tijdiUj, is transported by the fiuid: 

dtuj + {u ■ V)uj - vV^uj = F (3) 

with F = Sijdifj. The correlation function of the vorticity forcing term is thus: 

{Fix,t)Fiy,s))^G{x^y)S{t^s) (4) 

with G = — V^C. In particular G'(O) = 26^,. The fact that the correlation function of the vorticity forcing 
is a gradiant will have physical consequences. Physically eq.(^ means that for smooth solutions any powers 
of the vorticity, and in particular the enstrophy J ^, is conserved in absence of viscocity and forcing. 



Since the inviscid limit is of course not under analytical control, we have to make a few hypothesis which 
encode Kraichnan's scenario of inverse and direct cascades. The hypothesis are the following: 
(i) the velocity correlation functions are assumed to be smooth at finite viscocity and correlations of the 
velocity (without derivatives but at points coinciding or not) exist in the inviscid limit; 
(ii) Galilean invariant correlation functions, and in particular the velocity structure functions which are 
correlations of differences of the velocity field, are stationary; 

(iii) last but not least, in agreement with Kraichnan's picture, we demand that energy dissipative anomalies 
(but not enstrophy dissipative anomalies) be absent. 

The two first hypothesis are standard in statistical approach to turbulence while the third one is special to 
two dimensions. It follows by demanding that the enstrophy density fi = ^ is finite in the inviscid limit since 
the mean enstrophy density times the viscosity is equal to the mean dissipation rate, viTt) = ^((Vu) • (Vu)). 

Velocity correlations. Let us look at the two point velocity correlation function {u{x) ■ w(0)). As is 
well known, it satisfies the following equation at finite viscocity: 

dtiuix) • uiO)) - iv^((Au^-)(a;) {Auf{x)) + 2u{{Vu)ix) ■ (V7.)(0)) = d{x) (5) 

Here and in the following we shall denote velocity differences by (Aw'^)(a;) = u'^(x) — u'^(O). Eq.(P) assumes 
translation invariance and uses the fact that the external force is Gaussian and white-noice in time. Thanks 
to the fluid incompressibility the pressure drops out from this equation. The strategy consists in taking 
various limits of eq.(p|) in various orders. Let us take first the limit a: ^ followed by the inviscid limit. 
In this limit the second term in eq.dq) vanishes due to the assumed smoothness of the correlation functions, 
hypothesis (i). Recall now the hypothesis (iii) concerning the absence of energy dissipation. It in particular 
means that: 

lim limiy((VM)(x) •(Vu)(0)) =0 (6) 

Therefore the third term in eq.(P) also vanishes. It implies that in the inviscid limit the mean energy increases 
with time according to: 

a,(y)..o = ^5(0)^6 (7) 

Thus (^)iy=o — e^ up to a constant, and e is indeed the energy injection rate. This is simply the obvious 
statement that in absence of energy dissipation, and/or in absence of friction or other processes by which 
the energy may escape, all energy injected into the system is transfered to the fluid. It is expected to be 
transfered to the mode with the smallest possible momentum, the so-called condensate 0. In particular 
eq. (R) shows that in absence of energy dissipation a stationary state can not be reached although structure 
functions may converge at large time. This is one important difference between 2d and 3d turbulence. 

Let us now assume that the two point structure function is stationary, ie. dt{{Au)'^{x)) — 0, hypothesis 
(ii). From eq.(|5|) one obtains in the inviscid limit: 

iv^((Au'^)(x) (Au)2(x))..o = 26 - d{x) (8) 

Integrating it using parity invariance gives: 

{{Au^){x) {Au)'{x)),^o = 2 {-ex'' - O'^ix)) (9) 

with Va;-0(a;) = C{x). Eq.(^ together with eq.([l6|) below fully determine the three point velocity correlation. 
Although simple to derive, this equation seems not to have appeared in the existing literature n. 

Eq.(g) in particular shows that the inverse energy cascade takes place only if there is no dissipation 
anomaly, and thus only if the non-galilean invariant velocity correlation functions do not reach a stationary 
state (in absence of friction) . Of course this is also a direct consequence of the physical fact that the energy 
condensates into the mode of smallest possible momemtum. As expected, eq.(|9|) yields to Kolmogorov's 
scaling at large scale, since there Q''{x) vanishes, and Kraichnan's scaling at small scale since (ex'' — 0*'(a;)) ~ 
r^ at short distance. But one could be a little more precise. 

^ We thank G. Falkovich and R. Kraichnan for informations on the literature on turbulence. 



Let us first consider the short distance behavior in which the Kraichnan's direct cascade takes place. 
This corresponds to scale much smaller than the injection length a; << Li. There, 

Assuming isotropy and parity invariance the three point functions will be linear combinations of terms 
proportional to x^x^ x*' or {S'^^ x'' + 5^^x^ + 6^^x^)x'^ . Among these two proportionality coefficients only one 
of them could be fixed using eq.(|9|) only. However, the other coefficient is fixed by using another exact result 
for correlation functions mixing the vorticity and the velocity which is described in the following, see eq.(|l6[). 
This then gives for a; — > 0: 

((AM*)(x)(AuJ)(a;)(Au'=)(a;)) ~ ^ {{5'^x^ + ^^x' + 5^'x^)x^ - 2x'x^x^) (10) 
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For the transverse and longitudinal correlations this becomes: 

((A«)f) = ((Az.),|(Az.)i)c.+|ir3 (11) 

The coefficient e^, will be shown to be equal to the mean enstrophy dissipation rate in the following. Thus, 
as expected the 3-point velocity functions, which only depend on the enstrophy injection rate, are universal 
in the direct cascade. Eq.(pl]) may be called the "+1/8 law" of the direct cascade. 

Consider now the large distance behavior in which the Kraichnan's inverse cascade takes place, ie. scale 
X » Li. There, 

((Au'=)(a;) {/\uf{x))^2lx*' 

There is two possible terms for the three point functions, {5'^^x^ + 5^^x^ + 5^'^x^) and {x"^ x^ x^) / x^ whose 
coefficients can not be completely fixed using only cq.(g). But again these will be fixed by looking at the 
vorticity correlation functions, see below eq.([l^). One then gets for x — > cx): 

((Au*)(x)(AuJ')(x)(Au'=)(2;)) ~ Us'^x^ + P'^x' + S'^'x^) (12) 

with e the mean energy injection rate. Of course this gives the "+3/2 law" for the longitudinal statistics in 
the inverse cascade: 

((A^)i|) = 3((Au)||(A^)i)^+|r (13) 

This law and Kolmorgorov's scaling as well as the existence of a condensate in which the energy is accumu- 
lating was experimentally verified in |6|. 

Vorticity correlations. We now establish a formula for a mixed correlation function involving the 
vorticity and the velocity. Assuming that the structure functions of the velocity reach a stationary state 
implies that correlations of the vorticity also become stationary. The stationarity condition for the two point 
vorticity functions, ie. dt{uj{x)uj{y)) = 0, implies: 

- lv^((Au'=)(x) (A^)2(a;)) + 2^(V^(x) • V^(0)) = G{x) (14) 

with {Auj){x) — uj{x) — u!{0). As for the velocity correlations, let us first take the limit of coincident point 
a; — > at finite viscocity. The first term in eq. (^4|) then vanishes by the hypothesis on the smoothness of the 
correlation functions at finite v. Taking then the inviscid limit leads to: 

limz.((Vc.)2(x)) = iG(0) = e^ (15) 

This is just the usual statement on enstrophy dissipative anomaly. It equals the enstrophy dissipation rate 
and the enstrophy injection rate. Let us now take the limits in the reversed order, the inviscid limit first. In 
that case the second term of eq.(Q) drops out and one gets: 

-iv^((A^'=)(a;)(Ac.)2(a;)),=o-G(x) 



Assuming isotropy and parity invariance this gives: 

{{Au'^)(x) (Aujf{x)),=o = 2 — (^)(r) = -2e,„x'= + O(x') (16) 

r dr 

Thus this correlation is universal in the (UV) direct cascade regime. Its behavior at large scale depends on 
the way the forcing decreases at infinity. However the fact that these correlations decrease faster than 0(l/r) 
at infinity is linked to the fact the vorticity forcing correlation is a gradiant. The ultraviolet behavior (119) 
was also described in H. As mentioned above, this equation may be used to fix the coefficients of the infrared 
and ultraviolet expansions of the three point velocity function left undertermined by eq.(0). (One should 
use the relation 3((Au^)(Aa;)^) = 49^((Au^)'^) where u^ is the velocity component in complex coordinates 
z — X + iy.) But eq.(Qq) alone would not have been enough to determine these asymptotic expansions. 

Influence of anisotropy and parity symmetry breaking. Anisotropy is irrelevent both in the 
ultraviolet and in the infrared. Indeed, let us model anisotropy by incorporating higher spin components 
in the forcing correlation function C*-' (x) , assuming that they are still regular at the origin and decrease at 
infinity. These components will be subdominant in eq. (p|) both in the ultraviolet (since the spin n component 
will behave as r") and in the infrared (since they also vanish at infinity). 

Suppose now that parity symmetry may be broken. Eq.(P) is still valid (since it only assumes translation 
invariance) except that 6'^(a::) is determined up to 0*^ — * 9 + e''^djE{r). This may change the ultraviolet 
scaling of the transverse velocity correlations but not the scaling of the longitudinal correlations ((Au)^) 
although the amplitude may be modified. 

Influence of friction. In physical systems the infrared energy cascade will terminate at the largest 
possible scale at which the energy will escape. This could be mimic by introducing a friction term in the 
Navier-Stokes equation which then becomes: 

dtu^ + {u ■ W)u^ - v^'^v? + - li^' = -V^p + /^' (17) 

T 

with T the friction relaxation time, r > 0. Friction brings another inviscid characteristic length Lj into the 
problem: Lj ~ r^/^e^'^. It increases as the friction is reduced and one may suppose Li « Lj. It is the 
length at which the energy is extracted. The friction term dominates over the advection term at scale bigger 
than Lf. So the direct cascade should take place at distances x << Li << Lf and the inverse cascade at 
distances Li « x « Lf. 

Under the same hypothesis as before, the mean energy density relaxes in the inviscid limit according to 
dt{\) + ^{u^) — e. It therefore reaches a stationary limit with {u^) — It. Stationarity of the two point 
structure function, ie. 9t((Au)^(a;)) = 0, then gives in the inviscid limit: 

iv^((A^'=)(:r) (A«)2(x)) + i((A^)2(x)) = 26-e(a;) (18) 

Similarly, the stationarity of the vorticity correlations gives an equation similar to eq.(|l4) but with an 
extra term representing the friction. As for the case without friction let us first take the limit x ^ at 
finite viscocity and then the inviscid limit. Let us denote by Iw = ^^(0) the enstrophy injection rate 
and by Cw = hm^^o ^((Va;)^(a;)) the enstrophy dissipation rate. One then gets: {uj'^) = t (e^ ~?w) with 
(w^) = \iin^^o{oj'^{x)). It simply means that the enstrophy density is equal to the difference of the enstrophy 
injection and the enstrophy dissipation rates times the friction relaxation time. In particular, if r is finite so 
is the enstrophy density. As a consequence the vorticity two point correlation function {oj{x)uj{0)) will stay 
finite since it is bounded by (w^). Taking the limit in the reversed order, first v ^ 0, yields to the inviscid 
stationary equation: 

- iv^((Au'=)(x) {Ac^rix)) + ^(a>(x)c.(0)) = G{x) (19) 

Let us look at small distances in which the friction should be irrelevent. Let Q = liinx^o{Lo{x)uj{0)) , 
which is expected to be equal to T{ew — eiu) although nothing prevent it to be different. Since Q < re^, the 
second term in the l.h.s. of eq.([l9|) cannot dominate and V^((AM'')(a;) (A(jj)^(x)) ~ const, as x -^ 0. This 
implies that the velocity three point function {{Au)^) scale as r^. In other words, Kraichnan's scaling of the 
three point function is robust to friction in the direct cascade although the amplitude may change. 



More precisely, suppose that fi is finite and non- vanishing. Then the scahng formula (nO) for the three 
point function still holds but with e^j replaced by (e^r — ri)/r. Recall that it is likely that (e^r — fi)/T is 
equal to the enstrophy dissipation rate e^,, meaning that in presence of friction one has simply to replace 
the injection rate by the dissipation rate in formula dlQ). Moreover, the finiteness of the vorticity two point 
function at coincident points also imphes that: {{Au)^{x)) ~ § r^ since V^((Au)^(a;)) — 2{lj{x)uj{0)). It 
means that the UV energy spectrum is 

E{k)xTik-^ (20) 

without any logarithmic correction. This scaling may be broken only if fJ = li'ni,j;^Q{u;{x)u;(0))i,^Q vanishes. 
It is worth specifying more precisely in which scale domain eq.(pQ) will be valid. At finite viscocity, there 
are two ultraviolet characteristic lengths, the usual dissipative length r^ ~ i^'^/^ e^^'' and another friction 
length If ~ i/^'"^ T^'"^ above which friction dominates over dissipation. Since r^ << If in the limit we are 
considering, v ^ and t fixed, eq. (pO|) will be valid for If « x « Li. 

Let us now consider distances larger than the injection length Li but smaller than the friction length Lf. 
Then, 2e — C{x) ~ 2e and positivity argument cannot be applied. However, unless miraculous cancelations 
between the two terms in the l.h.s. of eq.(|8|), (which would mean that the domains in which advection 
or friction dominate intertwine), the correlation {{Au'^){x) (A'u)^(a;)) will still scale as r. But clearly this 
argument is less robust than the one used for the short distance analysis. 

Conclusions. Besides giving the expected formula for the 3-point velocity correlation functions, this 
short proof also indicates that if the inverse cascade takes place, as experimentally verified, then, in absence 
of friction, the non-Galilean invariant velocity correlation functions do not become stationary, although 
structure functions do. However, it gives no hints on how to decipher the behavior of the vorticity, one of 
the main challanging problems of two dimensional turbulence. 

Acknowledgements: We thank G. Falkovich and R. Kraichnan for informative correspondences. 
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